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Abstract. We review a three-dimensional formalism that provides a systematic way 
' to include relativistic effects including relativistic kinematics, the effects of negative- 

Q , energy states, and the boosts of the two-body system in calculations of two-body 

• bound-states. We then explain how to construct a conserved current within this rela- 

tivistic three-dimensional approach. This general theoretical framework is specifically 
CO , applied to electron-deuteron scattering both in impulse approximation and when the 

pTvy meson- exchange current is included. The experimentally-measured quantities A, 
B, and T20 are calculated over the kinematic range that is probed in Jefferson Lab 
experiments. The role of both negative-energy states and meson retardation appears 
to be small in the region of interest. 
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1 Introduction 



OO 

A number of the experiments being performed at the Thomas Jefferson Na- 
tional Accelerator Facility (TJNAF) involve the elastic and inelastic scattering 
of electrons off the deuteron at space-like momentum transfers of the order of 
q ' the nucleon mass. In building theoretical models of these processes, relativistic 

^ , kinematics and dynamics would seem to be called for. Much theoretical effort 

has been spent constructing relativistic formalisms for the two-nucleon bound 
state that are based on an effective quantum field theory lagrangian. If the 
usual hadronic degrees of freedom appear in the lagrangian then this strategy 
is essentially a logical extension of the standard nonrelativistic treatment of the 



■ two-nucleon system. 



Furthermore, regardless of the momentum transfer involved, it is crucial that 
a description of the deuteron be used which incorporates the consequences of 
electromagnetic gauge invariance. Minimally this means that the electromagnetic 
current constructed for the deuteron must be conserved. 

Of course, the two-nucleon bound state can be calculated and a correspond- 
ing conserved deuteron current constructed using non-relativistic NN potentials 
which are fit to the NN scattering data. This approach has met with consider- 
able success. (For some examples of this program see Refs. [1, 2].) Our goal here 
is to imitate such calculations — and, we hope, their success! — in a relativistic 
framework. To do this we construct an NN interaction, place it in a relativis- 
tic scattering equation, and then fit the parameters of our interaction to the 
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NN scattering data. We then calculate the electromagnetic form factors of the 
deuteron predicted by this NN model. By proceeding in this way we hope to 
gain understanding of the deuteron electromagnetic form factors in a model in 
which relativistic effects, such as relativistic kinematics, negative-energy states, 
boost effects, and relativistic pieces of the electromagnetic current, are explicitly 
included at all stages of the calculation. 

This program could be pursued using a four-dimensional formalism based on 
the Bcthc-Salpeter equation. Indeed, pioneering calculations of clcctron-deuteron 
scattering using Bethe-Salpeter amplitudes were performed by Zuilhof and Tjon 
almost twenty years ago [3, 4]. However, despite increases in computer power 
since this early work the four-dimensional problem is still a difficult one to solve. 
Since the NN interaction is somewhat phenomenological ultimately it is not clear 
that one gains greatly in either dynamics or understanding by treating the prob- 
lem four-dimcnsionally. Therefore, instead we will employ a three-dimensional 
formalism that incorporates what we believe are the important dynamical effects 
due to relativity at the momentum transfers of interest. 

We will use a three-dimensional (3D) formalism that, in principle, is equiva- 
lent to the four-dimensional Bethc-Salpcter formalism. This approach has been 
developed and applied in Refs. [5, 6, 7]. In this paper we will focus on the 
calculation of clastic clcctron-deuteron scattering. Here we review the formal- 
ism for relativistic bound states and show how to construct the corresponding 
electromagnetic current. Calculations of elastic electron-deuteron scattering are 
performed both in the impulse approximation and with some meson-exchange 
currents included. The results for the observables A, B and T 2 o are presented. 

Many other 3D relativistic treatments of the deuteron dynamics which are 
similar in spirit to that pursued here exist (see for instance Refs. [8, 9]). Of 
these, our work is closest to that of Hummel and Tjon [10, 11, 12]. However, 
in that work approximations were employed for ingredients of the analysis, 
such as the use of wave functions based on the 3D quasipotential propagator 
of Blankenbecler-Sugar [13] and Logunov-Tavkhelidze [14], approximate boost 
operators, and an electromagnetic current which only approximately satisfies 
current conservation. Calculations of elastic electron-deuteron scattering also 
were performed by Devine and Wallace using a similar approach to that pursued 
here [15]. Here we extend these previous analyses by use of our systematic 3D 
formalism. In this way we can incorporate retardations into the interaction and 
also use a deuteron electromagnetic current that is specifically constructed to 
maintain the Ward-Takahashi identites. 

The paper is organized as follows. In Section 2 we explain our reduction from 
four to three dimensions. In Section 3 we present a four-dimensional equation 
which is a modified version of the ladder Bethe-Salpeter equation. This mod- 
ified equation has the virtue that it, unlike the ladder BSE, incorporates the 
correct one-body limit. By applying our three-dimensional reduction technique 
to this four-dimensional equation we produce an equation which has the correct 
one-body limit and contains the correct physics of negative-energy states. In Sec- 
tion 4 we explain the various potentials that are used in calculations of deuteron 
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wave functions. These can be divided into two classes: instant potentials, and 
potentials that include meson retardation. Within either of these classes ver- 
sions of the potentials are constructed that do and do not include the effects 
of negative-energy states, in order to display the role played by such compo- 
nents of the deuteron wave function. Section 5 discusses our 3D reduction of the 
electromagnetic current that maintains current conservation. This completes the 
laying out of a consistent formalism that includes the effects of relativity system- 
atically, has the correct one-body limits, and maintains current conservation. In 
Section 6 we apply this machinery to the calculation of electron-deuteron scat- 
tering both in the impulse approximation and when corrections due to some 
meson-exchange currents are included. Finally, discussion and conclusions are 
presented in Section 7. 



2 The reduction to three dimensions 

The Bcthc-Salpctcr equation, 

T = K + KG T, (1) 

for the four-dimensional NN amplitude T provides a theoretical description 
of the deuteron which incorporates relativity. Here K is the Bethe-Salpeter 
kernel, and Go is the free two-nucleon propagator. In a strict quantum-ficld- 
theory treatment, the kernel K includes the infinite set of two-particle irreducible 
NN — > NN Feynman graphs. 

For the two-nucleon system an application of the full effective quantum field 
theory of nucleons and mesons is impractical and perhaps, since hadronic de- 
grees of freedom are not fundamental, inappropriate. In other words, the Bethe- 
Salpeter formalism may serve as a theoretical framework within which some 
relativistic effective interaction may be developed. But, if the NN interaction is 
only an effective one, then it would seem to be equally appropriate to develop 
the relativistic effective interaction within an equivalent three-dimensional for- 
malism which is obtained from the four-dimensional Bethe-Salpeter formalism 
via some systematic reduction technique. 

One straightforward way to reduce the Bethe-Salpeter equation to three di- 
mensions is to approximate the kernel K by an instantaneous interaction i^i ns t- 
For example, if q = (go,q) is the relative four-momentum of the two nucleons 
then 

K{q) = - r ^ - #(<!) = — 2T-2- ( 2 ) 

This, admittedly uncontrolled, approximation, yields from the Bcthc-Salpctcr 
equation the Salpeter equation: 

Tlnst = Ki nst + K[ nst (Go)Ti ns t , (3) 
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where the three-dimensional Salpeter propagator (Go) is obtained by integrating 
over the time-component of relative momentum, 



(Go> = J ^Go( P ;P). (4) 

Throughout this paper we denote the integration over the zeroth component of 
relative momenta, which is equivalent to consideration of an equal-time Green's 
function, by angled brackets. We shall consider only spin-half particles, and so 

(r \ - A t A 2 A 1 A 2 . m 

{Ga) E-e 1 -e 2 E + e 1 + e 2 > (5) 

where A ± are related to projection operators onto positive and negative-energy 
states of the Dirac equation, E is the total energy, and a = (pf + mf ) 1 / 2 . Note 
that for spin-half particles, this propagator (Go) is not invertible. 

In order to systematize this kind of 3D reduction one must split the 4D kernel 
K into two parts. One of these, K\ , is to be understood as a three-dimensional 
interaction in the sense that it does not depend on the zeroth component of 
relative four momentum 1 . We then seek to choose this K\ such that the 3D 
amplitude T\ defined by 

Ti = K x + Ki{G )T u (6) 

has the property that 

(Go)7i(G ) = (GoTGo). (7) 

It is straightforward to demonstrate that such a K i is denned by the coupled 
equations: 

K 1 = (G )- 1 (G Kg)(G )- 1 , (8) 
which is three-dimensional, and 

g = G Q + G Q {K-K 1 )g, (9) 
which is four dimensional. The K\ of Eq. (8) does this by ensuring that 

(G) = (Go). (10) 

The formalism is systematic in the sense that, given a perturbative expansion 
for the 4D kernel, K, a perturbative expansion for the 3D kernel, K\, can be 
developed. At second order in the coupling this gives: 

K[ 2) = (Go)- 1 (G ^ (2) Go)(Go)- 1 . (11) 

In H — | > H — h states this is just the usual energy-dependent one-particle-exchange 

interaction of time-ordered perturbation theory, but with relativistic kinematics, 
i.e. ignoring spin and isospin: 

~2 r 1 1 

(12) 



1 2uj 



1 +(!<-> 2) 



_E+ -ei-e' 2 -u 

1 Of course, this is not a covariant reduction, but covariance can be maintained by a 
suitable generalization of this idea [6] . 
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where u> is the on-shell energy of the exchanged particle. Note that (Go) must 
be invertible in order for the 3D reduction to be consistent. (Similar connections 
between three and four-dimensional approaches are discussed in Refs. [14, 16, 
17,18,19].) 

Equation (6) leads to an equation for the bound-state vertex function: 

A = K 1 (G )r 1 , (13) 

where i~i is the vertex function in the three-dimensional theory. The 4D vertex 
function, r, and the corresponding 3D one, A, are related via 

G r = gr 1 . (u) 



3 The one-body limit 

As mentioned above, and discussed many years ago by Klein [16], the propagator 
(Go) is not invertible and therefore the above reduction is not consistent. We 
shall show in this section that this difficulty is connected to the behavior of 
the three-dimensional equation in the one-body limit. In this limit we allow one 
particle's mass to tend to infinity. We expect that the amplitude T\ then reduces 
to that given by the Dirac equation for a light particle moving in the static field 
of the heavy particle. In fact, this does not happen unless we include an infinite 
number of graphs in the kernel of the integral equation Eq. (3). 

In fact, if a scattering equation with a kernel which contains only a finite 
number of graphs is to possess the correct one-body limit, two distinct criteria 
must be satisfied. First the 3D propagator should limit to the one-body propa- 
gator for one particle (the Dirac propagator in this case) as the other particle's 
mass tends to infinity. Second, as either particle's mass tends to infinity, the 
equation should become equivalent to one in which the interaction, K\, is static. 

Equation (3)'s lack of either of these properties stems from Eq. (1) not having 
the correct one-body limit if any kernel which does not include the infinite set 
of crossed-ladder graphs is chosen [20]. Solution of Eq. (1) with such a kernel is 
impractical in the NN system. Nevertheless, the contributions of crossed-ladder 
graphs to the kernel may be included in an integral equation for T by using a 
4D integral equation for K, the kernel of Eq. (1) 

K = U + UG C K. (15) 

Once Gc is defined this equation defines a reduced kernel U in terms of the 
original kernel K . The propagator Gc is chosen so as to separate the parts of 
the kernel K that are necessary to obtain the one-body limit from the parts 
that are not. U may then be truncated at any desired order without losing the 
one-body limits. The following 4D equation for the t-matrix is thus equivalent 
to Eqs. (1) and (15), 

T = U + U (G + Gc)T. (16) 
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We can now remedy the defects of our previous 3D reduction. Applying the 
same 3D reduction used above to Eq. (16) gives: 

Ti = (7i + Ui (Go + G C )T U (17) 

where the 3D propagator is 

A+A+ A+A 2 



{G0 + GC) - P° + -e 1 -e 2 -24-P0+ + e 1+ e 2 
A^A+ A-A- 



(18) 



and parameter that enters through the construction of Gc- This three- 

dimensional propagator was derived by Mandelzweig and Wallace with the choice 
k° 2 = P°/2 - (ml - ml)/(2P°) [21, 22]. With k\ chosen in this way (G + G c ) 
has the correct one-body limits as either particle's mass tends to infinity and 
has an invertible form. The kernel U\ is defined by Eqs. (8) and (10) with the 
replacements Go — * Go + Gc K — > U, and K\ — > U\. 

Here we are interested in the scattering of particles of equal mass and so we 
make a different choice for n 2 - Specifically, 

4 . £^li±S. (1 „ 

This form avoids the appearance of unphysical singularities when electron-deuteron 
scattering is calculated [7]. It yields a two-body propagator: 

{G ° + = P« + -e,-e 2 ----- po- +ei+£2 ■ ( 2 °) 

which is consistent with that required by low-energy theorems for Dirac particles 
in scalar and vector fields [23] . Another way of saying this is to realize that if we 
compare the the +H ► ++ piece of the amplitude 

^(Gq + Gc)^ (21) 

to the amplitude obtained at fourth order in the full 4D field theory then the 
contribution of negative-energy states agrees at leading order in 1/M [7]. 

For bound states the argument of the previous section leads to the 3D equa- 
tion: 

A = U 1 (Go + G c )r 1 . (22) 

Equation (22) is a bound-state equation which incorporates relativistic effects 
and the physics of negative-energy states. For instance, fig. 1 is one example of 
a graph which is included if Eq. (22), even if only the lowest-order kernel u[ 2 ^ 
is used, because of our careful treatment of the one-body limit. 
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Fig. 1. One example of a Z-graph which is included in our 3D equation (22). 

4 Results for the deuteron 

To calculate observables in the deuteron we now consider two types of kernels 
Ui , both of which are calculated within the framework of a one-boson exchange 
model for the NN interaction: 

1. Ui = J/instj the instantaneous interaction. 

(2) 

2. A kernel which is a retarded interaction. This is obtained from Eq. (11) 
by the substitutions K[ 2 ^ — > u[ 2 ^ and Go — > Go + Gc- 

These interactions are used in a two-body equation with the full ET Green's 
function given by Eq. (22), and also in an equation in which only the ++ sector 
is retained. For the instant interaction, we follow the practice of Devine and 

Wallace [15] and switch off couplings between the ++ and sectors, and 

between the H — and — + sectors. A partial justification of this rule follows from 
an analysis of the static limit of our 3D retarded interaction. 

The mesons in our one-boson exchange model are the 7r(138), the cr(550), 
the 77(549), the p(769), the w(782), and the (5(983). All the parameters of the 
model, except for the a coupling, are taken directly from the Bonn-B fit to the 
NN phase shifts [24] — which is a fit performed using a relativistic wave equation 
and relativistic propagators for the mesons. The a coupling is varied so as to 
achieve the correct deuteron binding energy for each interaction considered. Of 
course, we should refit the parameters of our NN interaction using our different 
scattering equations. However, for a first estimate of the importance of negative- 
energy states and retardation we adopt this simpler approach to constructing 
the interaction. Work on improving the NN interaction model is in progress [25]. 

Once a particular interaction is chosen, the integral equation (22) is solved 
for the bound-state energy. In each calculation, the a coupling is adjusted to get 
the correct deuteron binding energy, producing the results (accurate to three 
significant figures) given in Table 1. The value given for the instant calculation 
with positive-energy states alone is that found in the original Bonn-B fit. In all 
other cases the a coupling must be adjusted to compensate for the inclusion of 
retardation, the effects of negative-energy states, etc. We believe that this ad- 
justment of the scalar coupling strength is sufficient to get a reasonable deuteron 
wave function. The static properties of this deuteron are very similar to those of 
a deuteron calculated with the usual Bonn-B interaction. 

With the bound-state wave function in the center-of-mass frame has been 
determined in this fashion, it is a simple matter to solve the integral equation 
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(22) in any other frame. We choose to calculate electron-deuteron scattering in 
the Brcit frame. The interaction is recalculated in the Breit frame for a given Q 2 , 
and then the integral equation is solved with this new interaction. Because the 
formalism we use for reducing the four-dimensional integral equation to three 
dimensions is not Lorentz invariant there is a violation of Lorentz invariance in 
this calculation. Estimations of the degree to which Lorentz invariance is violated 
are displayed in Ref . [7] . 



Table 1. Sigma coupling required to produce the correct deuteron binding energy in 
the four different models under consideration here. 



Interaction 


States included 




Instant 


++ 


8.08 


Retarded 


++ 


8.39 


Instant 


All 


8.55 


Retarded 


All 


8.44 



5 Current conservation 

5.1 Currents in the three-dimensional formalism 

As discussed in the Introduction, we now want to compare the predictions of this 
formalism with experimental data gained in electron scattering experiments. In 
calculating the interaction of the electron with the hadronic bound state it is 
crucial to derive a 3D reduction of the electromagnetic current which is consistent 
with the reduction of the scattering equation we have chosen to use here. 

The current in the full four-dimensional formalism is obtained by coupling 
photons everywhere on the right-hand side of Eq. (1). This produces the following 
gauge-invariant result for the photon's interaction with the bound state: 

A„ = r(P')Go(P')J^G (P)r(P) 

+ r(P')G Q (P')K~<G Q (P)r(P), (23) 

where P and P' are the initial and final total four-momenta of the deuteron 
bound state. Here contains the usual one-body currents and represents 
two-body contributions which are necessary for maintaining the Ward-Takahashi 
identities. All integrals implicitly are four-dimensional. The connection to the 
three-dimensional amplitude, P\, obtained from Eq. (22) is made by inserting 
Eq. (14) into Eq. (23), giving 



a, = r^P'mp') [J, + k;] s(p)>A(p). 



(24) 
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Once the effective operator {G(P') [Jp + K^j Q{P)) is calculated the expression 
(24) involves only three-dimensional integrals. 

Since Q is an infinite series in K — K\ this result would not be much help 
on its own. But, given a result for F\ obtained by systematic expansion of Ki, 
the amplitude A^ can be analogously expanded in a way that maintains current 
conservation. K\ as defined by Eq. (8) is an infinite series and the condition 
(10) is imposed order-by-order in the expansion in K — K\ defines K\ to some 
finite order. The question is: Does a corresponding 3D approximation for the 
current matrix element (24) exist that maintains the Ward-Takahashi identities 
of the theory? It turns out that the current matrix element (24) is conserved if 
G{Jn + K^)Q on the right-hand side of Eq. (24) is expanded to a given order 
in the coupling constant and the kernel K\ used to define i~i is obtained from 
Eq. (10) by truncation at the same order in the coupling constant. 

This is done by splitting the right-hand side of Eq. (24) into two pieces, one 
due to the one-body current J M , and one due to the two-body current K^. If Ki 

has been truncated at lowest order — i.e., K\ — k[ 2 ^ — then, in the J M piece, we 
expand the Qs and retain terms up to the same order in — k[ 2 \ A piece 
from the two-body current, in which we write Q = Go, is added to this. That is, 

(2) 

we define our second-order approximation to A^, A^ , by 

4 2 ) = r 1 (P , )<G^)A(P) 

+ r 1 (P')(G (P')(K^(P')-K[ 2 \p'))G^)r 1 (P) 
+ r 1 (P'){G% ll (KW(P) ^ 2) (P))Go(P))A(P) 
+ r 1 {P , ){Go(P')^Go{P))r 1 {P), (25) 
where Gq^ = Gq(P')J IJi Go(P). It can now be shown that if Eq. (10) expanded 

(2) 

to second order defines K\ , the corresponding amplitude for electromagnetic 
interactions of the bound state, as defined by Eq. (25), exactly obeys 

Q»A {2) = 0. (26) 

It is straightforward to check that the same result holds if Eq. (10) for K\ is 
truncated at fourth order, while the one-body and two-body current pieces are 
expanded to fourth order. 

(2) 

The amplitude ^ includes contributions from diagrams where the photon 
couples to particles one and two while exchanged quanta are "in-flight" . These 
contributions are of two kinds. Firstly, if the four-dimensional kernel K is depen- 
dent on the total momentum, or if it involves the exchange of charged particles, 
then the WTIs in the 4D theory require that contain terms involving the 
coupling of the photon to internal lines in K. Secondly, even if such terms are 
not present, terms arise in the three-dimensional formalism where the photon 
couples to particles one and two while an exchanged meson is "in-flight" . These 
must be included if our 3D approach is to lead to a conserved current. (See Fig. 2 
for one such mechanism.) 
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Fig. 2. One example ol a two-body current that is required in our formalism in order 
to maintain current conservation. 

A special case of the above results occurs when retardation effects are omit- 
ted, i.e., the kernel K\ — Ki DSt , is chosen, and the bound-state equation (13) is 
solved to get the vertex function F\ = il ns t- Then a simple conserved current is 
found: 

Anst^ = r inst (p')(<^>r inst (p) + r inst (p')(Go(P'))Kl st ^G a (p))r inst (p) 7 

(27) 

where we have also replaced the meson-exchange current kernel by the instant 
approximation to it. 

5.2 Current conservation in the 4D formalism with Gc 

In Ref. [7] we showed how to construct a conserved current consistent with the 
4D equation 

r = U(G + G c )r. (28) 

This turns out to be a moderately complicated exercise, because the propagator 
Gc depends on the three-momenta of particles one and two, not only in the 
usual way, but also through the choice (19) made for above. However, a 4D 
current Qq ^ = Gq^ + corresponding to the free Green's function Go + Gc 
can be constructed. Its form is displayed in Ref. [7] and is not really germane 
to our purposes here, for, as we shall see hereafter, only certain pieces of the 
current Qq ^ are actually used in our calculations. 

5.3 Reduction to 3D and the ET current 

Having constructed a 4D current for the formalism involving Gc that obeys 
the required Ward-Takahashi identity, we can apply the reduction formalism of 
Section 5.1 to obtain the currents corresponding to the 3D reduction of this 4D 
theory. The result is: 

4 2) =A,et(P')(^)A,et(P) 

+ A,et(p')((Go + Gc)(p')(k^{p') u[ 2 \p'))g^)r hET (P) 
+ r hET (P')(Gl,(K^(P) u[ 2 \p))(g + G c )(P))r hET (P) 

+ A,et(P')((Go + G c )(P')k;W(Go + G c )(P))A,et(P), (29) 
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where A,et is the solution of Eq. (22) with U\ — U\ ' . This current obeys 
the appropriate Ward-Takahashi identity. In fact in one-boson exchange models 
the only contributions to K^ 2 ' give rise to isovector structures, and so their 
contribution to electromagnetic scattering off the deuteron is zero. 



5.4 Impulse-approximation current based on the instant 
approximation to ET formalism 

Just as in the case of the Bethe-Salpeter equation, if the instant approximation is 
used to obtain a bound-state equation with an instant interaction from Eq. (28) 
then a corresponding simple conserved impulse current can be constructed: 

"4inst,/i = ^lnst(^0,/i)-^inst- (30) 

Now we note that the full result for Gq „ was constructed in order to obey 
Ward-Takahashi identities in the full four-dimensional theory It is not necessary 
to use this result if we are only concerned with maintaining WTIs at the three- 
dimensional level in the instant approximation. Therefore we may construct the 
corresponding current 

G?n S t,M^;P, Q) = i(d 1 (p 1 )d 2 ( P 2 + Q)3 {2) d 2 { P 2) 

+d 1 ( Pl )4( P2 + Q)jWd c 2 (p 2 )) + (1 <- 2). (31) 

Here di is the Dirac propagator for particle i, and — q"/^ is the usual one-body 
current, with q is the charge of the particle in question. Meanwhile d^ is a one- 
body Dirac propagator used in Gc(P) to construct the approximation to the 
crosscd-laddcr graphs. Correspondingly d\ appears in Gc(P + Q), which does 
not equal d\, even if particle i is not the nucleon struck by the photon. Finally 

;$ = ?27,-if, (32) 

where 

- ( 2) ?V +P2 P . . 

J l = 92— r— — 72o, (33) 
e 2 + £ 2 

with p 2 = (e(p2),P2)- (For further explanation of these quantities and the ne- 
cessity of their appearance here the reader is referred to Ref. [7].) 

If a vertex function J] nst is constructed to be a solution to Eq. (22) with an 
instant interaction then the three-dimensional hadronic current: 

-^inst./j -Anst^i ns t^-^inst (34) 

is conserved. This current is simpler than the full ET current and omits only 
effects stemming from retardation in the current. Our present calculations are 
designed to provide an assessment of the role of negative-energy states and re- 
tardation effects in the vertex functions. Therefore we use the simple current 
(34) in all of our calculations here — even the ones where A is calculated using 
a retarded two-body interaction. The effects stemming from retardation in the 
current are expected to be minor, and so we expect this to be a good approx- 
imation to the full current in the three-dimensional theory. Future calculations 
should be performed to check the role of meson retardation in that current. 
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6 Results for electron-deuteron scattering 
6.1 Impulse approximation 

We arc now ready to calculate the experimentally observed deuteron electro- 
magnetic form factors A and B, and the tensor polarization T20. These are 
straightforwardly related to the charge, quadrupolc, and magnetic form factors 
of the deuteron, Fc, Fq, and Fm- These form factors in turn are related to the 
Brcit frame matrix elements of the current *4 M discussed in the previous section, 

F C = 1 (<0|^°|0) + 2{+l\A {) \ + 1», (35) 

F Q = — -L_ ( <0 |.A°|0) - (+1|^| + 1)), (36) 

Fm = ~\ (+11-4+10), (37) 

where | + 1), |0) and | — 1) are the three different spin states of the deuteron. 

We take the wave functions constructed for the four different interactions of 
Section 4 and insert them into the expression (34). In using any of the interac- 
tions obtained with only positive-energy state propagation we drop all pieces of 
the operator G? nst „ in negative-energy sectors. 

The single-nucleon current used in these calculations is the usual one for 
extended nucleons. We choose to parametrize the single-nucleon form factors 
F\ and F2 via the 1976 Hohlcr fits [26]. Choosing different single-nucleon form 
factors does not affect our qualitative conclusions, although it has some impact 
on our quantitative results for A, B, and T20. 

Using this one-body current we then calculate the current matrix elements 
via Eq. (24). This is a conserved current if the vertex function ri is calculated 
from an instant potential. However, if a potential including meson retardation 
is used it violates the Ward-Takahashi identities by omission of pieces that are 
required because of the inclusion of retardation effects in the calculation. Work 
is in progress to estimate the size of these effects. 

The results for the impulse approximation calculation of the experimental 
observables A, B, and T20 are displayed in Fig. 3. We also show experimental 
data from Rcfs. [27, 28, 29, 30, 31] for A, from Refs. [29, 30, 32, 33] for B and 
from Ref. [34] for T20. A number of two-body effects must be added to our cal- 
culations before they can be reliably compared to experimental data. However, 
even here we see the close similarity of the results for these observables in all 
four calculations. The only really noticeable difference occurs at the minimum 
in B. There, including the negative-energy states in the calculation shifts the 
minimum to somewhat larger Q 2 . A similar effect was observed by van Ordcn 
et al. [8] in calculations of electron-deuteron scattering using the spectator for- 
malism. However, note that here, in contradistinction to the results of Ref. [8], 
the inclusion of negative-energy states does not bring the impulse approximation 
calculation into agreement with the data. 
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Fig. 3. The form factors A(Q 2 ) and B(Q 2 ) and the tensor polarization T20 for the 
deuteron calculated in impulse approximation. The dash-dotted line represents a cal- 
culation using a vertex function generated using the instant interaction. Meanwhile 
the solid line is the result obtained with the retarded vertex function. The dotted and 
long dashed lines are obtained by performing a calculations with instant and retarded 
interactions in which no negative-energy states are included. 



The fact that negative-energy states seem to have a smaller effect on ob- 
servables in the ET analysis than in the spectator analysis of van Ordcn ct 
al. [8] is somewhat surprising since our "ET" propagator has twice the negative- 
energy state propagation amplitude of the spectator propagator. Thus, other 
differences between the ET and spectator models, not just differences in the role 
of negative-energy states in the two approaches, appear to be responsible for 
Ref. [8]'s success in reproducing the minimum in B. 

For the tensor polarization T20 the different models produce results which are 
very similar. This suggests that this observable is fairly insensitive to dynamical 
details of the deuteron model, at least up to Q 2 = 4 GeV 2 . 



6.2 Meson-exchange currents 

As Q 2 increases the cross-section due to the impulse approximation diagrams 
drops precipitously. Thus we expect that in some regime other interactions may 
become competitive with the impulse mechanism. One such possibility is that 
the photon will couple to a meson while that meson is in flight. Because of the 
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deuteron's isoscalar nature and the conservation of G-parity, the lowest mass 
state which can contribute in such meson-exchange current (MEC) diagrams is 
one where the photon induces a transition from a n to a p. 

This pirj MEC is a conserved current whose structure can be found in 
Rcfs. [10, 35]. The couplings and form factors for the meson- nucleon-nucleon 
vertices are all taken to be consistent with those used in our one-boson-exchange 
interaction. Meanwhile, the pnj coupling is set to the value g p7I1 — 0.56, and a 
vector meson dominance form factor is employed at the pirj vertex: F p ^ 1 {q) = 
l/(q 2 — m^,)- The value of this MEC is added to the impulse contribution cal- 
culated above and A, B, and T 2 o are calculated. This is done with the vertex 
function obtained from an instant interaction, and consequently the electromag- 
netic current is exactly conserved. The results of this calculation are displayed 
in Fig. 4. We see that at Q 2 of order 2 GeV 2 the pn-y MEC makes a significant 
contribution to all three observables. However, far from improving the agreement 
of the position of the minimum in the B form factor with the experimental data, 
this particular MEC moves the theoretical result away from the data — as noted 
by Hummel and Tjon [10], and seen within a simplified version of the formalism 
presented here by Devine [35]. Thus, it would seem that some physics beyond 
the impulse approximation other than the pirj MEC plays a significant role in 
determining the position of the minimum in B(Q 2 ). 




Fig. 4. The form factors A(Q 2 ) and B(Q 2 ) together with the tensor polarization for 
the deuteron. The long dashed line is an impulse approximation calculation with an 
instant interaction. The solid line includes the effect of the pny MEC. 
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7 Conclusion 

A systematic theory of the electromagnetic interactions of relativistic bound 
states is available in three dimensions. In this formalism integrations are per- 
formed over the zeroth component of the relative momentum of the two particles, 
leading to the construction of "equal-time" (ET) Green's functions. If the for- 
malism is to incorporate the Z-graphs that are expected in a quantum field 
theory, then the propagator must include terms coming from crossed Feynman 
graphs. Here we have displayed a three-dimensional propagator that includes 
these effects correctly to leading order in 1/M. 

Given a suitable choice for the ET propagator, the electromagnetic and in- 
teraction currents which should be used with it can be calculated. If these are 
truncated in a fashion consistent with the truncation of the NN interaction in 
the hadronic field theory then the Ward-Takahashi identities are maintained in 
the three-dimensional theory. A full accounting of the dynamical role played by 
negative-energy states and of retardations in electromagnetic interactions of the 
deuteron is thereby obtained. 

Calculations have been performed for both the impulse approximation and 
when the p-wy MEC is included. In our MEC calculations we use an instant 
approximation for the electromagnetic current. This current satisfies current 
conservation when used with deuteron vertex functions that are calculated with 
instant interactions. We also have used this simpler current with vertex func- 
tions which are calculated with the retarded interactions obtained within the 
ET formalism. 

Comparing impulse approximation calculations with and without negative- 
energy states indicates that the role played by negative-energy state components 
of the deuteron vertex function is small. This corroborates the results of Hum- 
mel and Tjon and is in contrast to those obtained in Ref. [8]. Because the ET 
formalism incorporates the relevant Z-graphs in a preferable way, we are confi- 
dent that these Z-graphs really do play only a minor role in calculations that 
are based upon standard boson-exchange models of the NN interaction. 

The results for impulse approximation calculations of the electromagnetic ob- 
servables are relatively insensitive to the distinction between a vertex calculated 
with retardations included and one calculated in the instantaneous approxima- 
tion. The results of both calculations fall systematically below experimental data 
for the form factors A and B for Q of order 1 GeV. This deficiency at higher Q 
suggests that mechanisms other than the impulse approximation graph should 
be significant. Indeed, when the pnj MEC graph is included in our calculation 
it somewhat remedies the result for A(Q 2 ). However, it fails to narrow the gap 
between our result for B(Q 2 ) and the existing experimental data. The signifi- 
cant gap that remains between our theoretical result for B(Q 2 ) and the data 
indicates that it is an interesting observable in which to look for physics of the 
deuteron other than the simple impulse mechanism or the standard pirj MEC. 
Finally, the existing tensor polarization data are reasonably well described. This 
is consistent with previous analyses which have shown T 2 o to be less sensitive to 
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non-impulse mechanisms. 
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